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We develop a scaling theory for the finite-size critical behavior of the microcanonical entropy 
(density of states) of a system with a critically-divergent heat capacity. The link between the 
microcanonical entropy and the canonical energy distribution is exploited to establish the former, 
and corroborate its predicted scaling form, in the case of the 3d Ising universality class. We show that 
the scaling behavior emerges clearly when one accounts for the effects of the negative background 
constant contribution to the canonical critical specific heat. We show that this same constant plays 
a significant role in determining the observed differences between the canonical and microcanonical 
specific heats of systems of finite size, in the critical region. 
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I. INTRODUCTION 



Statistical mechanics can be formulated in any of a set of ensembles distinguished by the relationship between the 
system and its environment ^ . The principal members of this set are the microcanonical (prescribed energy) and 
canonical (prescribed temperature) ensembles. In the thermodynamic limit (when it exists) the ensembles yield the 
same predictions (and are, in this sense, equivalent) and the choice of ensemble is a matter of practical convenience. 
The canonical ensemble tends to win this contest because it circumnavigates the hard-constant-energy constraint 
imposed by the microcanonical ensemble. 

The two ensembles are, however, not always equivalent They differ for systems which are 'small' in some 
sense: inherently small systems such as nuclei or clusters Q ; systems with unscreened long-range forces where the 
thermodynamic limit is problematic; and systems at critical points which are our principal concern here. 

Theoretical studies of critical phenomena are almost invariably conducted within the framework of the canonical 
ensemble In consequence there is no substantive framework within which to interpret computational studies 
of microcanonical critical behavior. Such studies do, nevertheless, exist, having been motivated, variously, by the 
belief that the microcanonical framework may have some computational advantages and by the discovery |^ that, 
apparently, critical anomalies in the microcanonical heat capacity are significantly enhanced with respect to their 
canonical counterparts. 

This paper goes some way towards supplying the missing framework. We develop (section 2) a finite-size-scaling 
theory |^ for the microcanonical entropy (the density of states) of a system with a critically-divergent heat capacity. 
In so doing we have, of necessity, to consider more general questions about the structure of the density of states of a 
finite-size system - in p articular the implications of well-established results for the finite-size structure of the canonical 
partition function flCt] . 




Though somewhat more than a phenomenology, our theory falls short of being microscopically explicit: to determine 
an explicit form for the relevant scaling function we need to appeal (section 3) to Monte Carlo (MC) measurements 
of the critical canonical energy probability distribution (pdf). 

The canonical energy pdf itself has a near-critical finite-size-scaling form which has featured in a number of studies 
of critical points in fiuids and lattice gauge theories . Since energy fluctuations (like the critical anomaly in 
the canonical specific heat which they control) are relatively weak (by comparison with the fluctuations of the order 
parameter, and the divergence of its response function) the degree of 'scaling' reported in previously measured energy 
pdfs has been relatively poor -unsatisfactorily so for our purposes here. This problem is addressed in section 3. We 
show that one can fold out (from the measured distributions) the sub-dominant (but signiflcant) non-scaling effects 
that are associated with the constant background contribution to the canonical heat capacity, negative in the case of 
the 3d Ising model This procedure exposes the underlying behavior, which manifests scaling to an impressive 
degree. In addition to providing us with the platform needed for this work, this procedure may offer the basis for 
improving the mixed-scaling-fleld theory [pT| of critical points in systems that belong to the Ising universality class 
but which do not have full Ising symmetry; recent studies have suggested that the current framework is not fully 
satisfactory. 

The scaling form for the critical energy pdf allows us to determine the scaling form of the microcanonical entropy. In 
section 4 we explore this form and show that it is consistent with predictions for both the bulk-critical limit (as regards 
the parameters characterizing the speciflc heat singularity p3[|) and the finite-size critical limit (the Fisher-Privman 
constant Q). 

The microcanonical entropy also provides us with a unified basis for dealing with both the canonical and the 
microcanonical specific heats (section 5). We show that the 'corrections' to the scaling behavior of the canonical specific 
heat (the negative background constant) have subtle consequences for the microcanonical behavior. In particular they 
serve to amplify the difference between microcanonical and canonical behavior, and are at least partially responsible 
for the strength of the anomaly observed in some microcanonical studies [|| . 



We consider a d-dimensional many-body system of linear dimension L; we assume hypercubic geometry with periodic 
boundary conditions. The canonical partition function is, in principle, a discrete sum over system microstates (r) or 
system energy levels (s): 




A. The microcanonical scaling ansatz 
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where fi^ is the degeneracy of level s. We shall suppose that the system is sufSciently large that the sum over levels 
can be replaced by an integral: 



Z(/?,L) = / deVL(e,L)e 



(lb) 



where e = E/L"^ is the energy density. The function fl{£, L) is the density of states; as we have defined it, it is a true 
density, having dimensions of inverse energy. We note that the transition from the discrete representation (Eq. [l^) to 
its continuum counterpart (Eq. lb) requires some care: it is discussed in Appendix A. 

Our microcanonical scaling theory comprises a proposal for the form of the density of states function. Wc formulate 
it in two stages. Consider first a regime remote from critical points or lines of phase coexistence. In such a regime we 
make the general finite-size ansatz [nSl: 
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The structure proposed for the prefactor makes this a little more than simply a definition of the microcanonical 
entropy density s(e). In its support we note, first, that one may readily verify it explicitly (Appendix B) in the case 
of some simple model systems. Secondly we note the implications for the associated canonical partition function. 
Inserting Eq. (^ into Eq. ( p^ , a saddle-point integration gives 



Z(/3,L) 





"7 


2tt_ 





where 



and e is the solution of 



de[-s"{e)]'^\^"^<^^-P^ 

fi(3) =/3e~ s{i) 
/3 = s'{e) 



(3) 
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Eq. (|^) recovers the prefactor-free form of the canonical partition function believed to be widely appropriate in regions 
(those where the saddle point integration is to be trusted) remote from critical points or lines of phase coexistence 
[Q. We note that this form is achieved by virtue of the prefactor that does appear in the density of states ansatz 
(Eq. which is just such as to cancel the contributions made by the fluctuations about the saddle point [ p^ . 

The argument we have given leaves open the possibility of power-law corrections to Eq. (^. It has long been 
believed, and more recently established rather generally |Q, that the corrections to the leading form are actually 
exponentially small in the system size. Since the saddle-point integration necessarily generates power-law corrections, 
one must suppose that there are compensating power-law corrections to the ansatz (Eq. |^) for the density of states. 
This conclusion serves as a warning (already suggested by the double appearance of the function s(e) in Eq. (^)) that 
the microcanonical framework faces problems which are skirted in the canonical formalism jl^ . 

Now, more specifically, consider a system, of the kind specified above, in the vicinity of a critical point. We will 
suppose that the critical point has a divergent heat capacity; where we need to be more specific we shall assume 
it is a member of the d=3 Ising universality class (or, more specially still, the d=:3 Ising model itself). Within the 
microcanonical framework the critical point of such a system is located by a critical value £c of the energy density, 
sharply-defined in the thermodynamic limit. We are concerned with the behavior of the microcanonical entropy for 
energies in the vicinity of this critical value. To describe this regime we introduce the dimensionless scaling variable 



X = a^L^^'^" (e — Ec) 
where Oe is an appropriate scale factor and the index is defined by 
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with a the index (assumed positive) characterizing the heat capacity divergence, and v the correlation length index 
p9| . We now reformulate and extend our basic ansatz (Eq. |^) with the proposal that, in a region of sufficiently large 
L and sufficiently small | e — Cc | 
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with 



L'^.3{e,L) 



+ /3c(e-ec)]+5(a;) 
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Here Sc is an unimportant constant, Pc is the critical inverse temperature and S{x) is a finite-size-scaUng function, 
universal given some convention on the scale factor a^, introduced in Eq. (^. The remainder of this paper is devoted 
to providing support for this proposal, and exploring the structure of the microcanonical entropy scaling function 
which it introduces. 



II. DETERMINING THE SCALING FUNCTION 

It should be possible to determine the finite-size-scaling function S{x) within the renormalization group framework 
We have not done that. Instead we have chosen to learn what we can about this function from its signatures in 
MC studies of the canonical ensemble. 

Consider, then, the implications of the scaling form Eq. (|8b|) for the canonical partition function, Eq. (p^). We 
suppose initially (we shall have to refine the supposition, shortly) that the relevant part of the energy spectrum is 
adequately captured by Eq. (|^). Then 

Z(/3,L)^e-^'/«('')i(y) (9) 

where 

/o(/3) =/3ec-Sc (10) 

and 



Z[y)^ J dJ ^^e'^y^^^^^ (11) 



while 

y = a:'L'/^f3-Pc) (12) 

provides a scaling measure of the deviation from the critical temperature. We have made use of the hyperscaling 
relation [Tg| ] which links the correlation length index v and the heat capacity index a through 

1 1 2-a , , , 

- + — = = d (13) 

The scaling form of the free energy follows: 

F{i5, L) = - \nZ{P, L) ~ - Ini(y) = F{fi, L) (14) 

The canonical energy pdf 

P(e|/3,L)^Z-i(/3,L)f7(e,L)e-''^'^ (15) 

may also be written in scaling form: 

P{e\(3,L)de = P{x\y,L)dx (16) 

with 



P{x\v, L) ^ Z~\y)J ^Me-^+'^C-) ^ p^^\y) (17) 
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The scahng predictions for the pdf may be tested by examining its cumulants P4 I, for which the free energy is a 
generator: 



(/3,L)^(-1) 
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Eq. (|lj) then impUes that the cumulants have the scahng form 



1 —Ti 



where the scaled cumulants x*^"^ (y) are universal functions: 

The canonical mean of the energy density at criticality (/? = /3c) follows as: 



(18) 



(19) 



(20) 



(21) 



MC measurements on the 3d Ising model using a range of system sizes (Fig. 1) are fully consistent with this behavior. 
Eq. ( [1^ ) implies, likewise, that the canonical variance of the energy density should have the power law behavior 



(22) 



MC measurements (Fig. 2) are only partially consistent with this prediction: the power law is confirmed, but with 
an extrapolation whose intercept is far from zero. This inconsistency is reflected in the rather limited success (Fig. 
3) of attempts to collapse the measured energy pdfs for different system sizes on to a single scaling form. The source 
of these problems can be guessed from the implications of Eq. ( p^ ) for the canonical specific heat, which it mirrors: 
the scaling form fails to capture the effects associated with the constant background which constitutes the dominant 
correction to pure scaling (power-law divergence) of the canonical specific heat. 

There are two ways to rectify this failure. One might extend the theory to predict the behavior of the (very) finite 
systems accessible to MC study; or one might seek to correct the MC results to expose the true scaling behavior. We 
adopt the latter strategy. 

Define 



A/(/3,L) = i-'' F{(3,L)-F{(3,L) 



(23) 



the difference between the true free energy density and its asymptotic scaling form (see Eq. (|14|)). We shall ignore the 
effects of confluent singularities: they are not the dominant 'corrections to scaling' here. Then Af{(3, L) is analytic 
and may be approximated, near f3c, by the expansion 



A/(/3,L)c.^A/, 



(24) 



n=0 



These additional contributions to the free energy imply additional contributions to the energy cumulants (Eq. |T^) : 

(25) 

At criticality Eq. ( |l9| ) must then be modified to read 



Ae(")(/3,L) ^ (_i)«+i^-(«-i)<i^!MM) 
^ ' ^ ' 9/3" 
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i(")(y = 0) + Ax(")(L) 
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where 



A4")(L) = x(")(y = 0,L) - i(")(2/ = 0) = (-l)"+i<'L'^-"/''A/i") 



(27) 
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The n = 1 correction is absent by fiat: the choice of Cc ensures this. The n > 3 corrections are sufficiently strongly 
'irrelevant' (they vanish sufficiently strongly with L) that they may reasonably be neglected. But the n = 2 correction 
decays only slowly: 



where the last step defines g{L) (a convenient parameter) while 
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is identifiable as the constant 'background' to the near-critical canonical specific heat. With this addition, Eq. (22) is 
modified to read 



[{e^), - {e)l] ^ L'^/'^a:^x(^Hy = 0, L) = L'^^'^aJ^xi'Hy ^ 0) + 



COc 



(30) 



which is now fully consistent with the MC measurements of Fig. 2, with (it is to be noted) a negative value for cqc [ p3[ . 
From a thermodynamic point of view these results simply reflect the fact that, for any system size practically accessible 
to computer-simulation, the 'critical' contribution to the canonical specific heat is not large enough to dominate the 
'non-critical background'. But the argument also shows us how to eliminate the effects of this 'background' from the 
energy pdf. 

Consider the cumulant representation |22] of the scaling energy pdf (Eq. n% at criticality: 



1 r°° 

P{x\y = 0) — / dre'^^ exp 
27r 



The corresponding relation for the observed energy pdf at criticality, written in its inverse form, is 



exp 
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(2/ = 0,L) 



dx'e-"^P(x'|j/ = 0,L) 
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Appealing to the our conclusion that, for large enough L, the cumulants of the two pdfs differ significantly only in 
the n = 2 case, and using Eqs. (p7| ) and (28), we find that 



2,71 



dr 



^J2^)J- 



dx'e'^''-'''^^-sW^/^P{x'\y = 0,L) 
(ia;'e-(^-^')'/[2''(^)lp(a;'|y = 0, L) 



(33) 



This result shows that the scaling form of the critical pdf may be exposed by convolution of the observed (and thus, 
generally, non-scaling) pdfs with gaussians whose widths are controlled by the specific heat background. Note that 
the argument rests on the fact that this background is negative (so that g{L) as defined in Eq. ( p8| ) is positive). If 
the background constant were positive our argument would have to be restructured to prescribe the scaling form by 
a process of deconvolution, which is numerically problematic. As it is, the convolution process can be implemented 
easily. With cqc fixed by the ordinate intercept in Fig. 2, the pdfs measured on different system sizes can each 
be corrected in this way to yield estimates of the scaling pdf. The results are shown in Fig. 4. The improvement 
with respect to the raw data (Fig. 3) is striking. This improvement reflects not only the removal of the non-scaling 
contribution to the second cumulant but also that the requisite convolution process provides a natural smoothing 
of the MC data ||2^. The consequences of this correction for the shape of the distribution are also striking. The 
skewness p9[ | clearly visible in the raw distributions (Fig. 3) is largely suppressed to expose a scaling form that is, at 
first appearance, gaussian. Indeed the portion of the distribution evident on the scale of Fig. 4 is gaussian to within 
deviations of a few per cent. However the behavior in the wings (evident on the logarithmic scale utilized in Fig. 5) 
is markedly different on the high- and low-energy sides. 

The scaling of the critical energy pdf corroborates the scaling of the microcanonical entropy (cf Eq. (p7|)). Given 
the double appearance of S{x) in Eq. ( [l7| ) it is practical to infer only the 'effective' microcanonical entropy scaling 
function 



Seffix) =S{x) + ^\n 



^S"{x) 
2n 



(34) 
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Fig. 5 shows the form imphed by Eq. (|r 

^eff{x) = Seffix = 0) +ln 



(35) 



P{x = 0\y = 0) 

We note as a matter of empirical fact that Seff{x) is concave. The concavity of S{x) itself is already assumed in our 



basic scaling ansatz |15 



III. THE SCALING THEORY: IMPLICATIONS AND TESTS 



Although we have no first-principles calculation of the scaling function S{x) to offer here we can identify, and test, 
some of the properties it must have to match anticipated behavior in both the thermodynamic and finite-size-critical 
limits. 

We consider, in particular, the limiting large | x \ behavior. In this regime we anticipate that 



S{x) ~ -6± I X I** +r± (I X |> 1) 



(36) 



where the -I- and — subscripts refer, respectively, to the regions of positive and negative x. To make explicit identi- 
fications of the new quantities introduced in this equation (the exponent 9 and the amplitudes h±, r±)^ consider the 
scaling part of the partition function (Eq. ^T]) . In the limit of large | y \ the integral in Eq. ( pi] ) is dominated by one 
or other of the large | x \ regimes. Substituting Eq. (|3^), a saddle-point integration yields 



lnZ(y):^a± lyr/C'-D+ri (| y |» 1) 
where the -I- and — subscripts now refer, respectively, to the regions of negative and positive y 



and 



a+ 



1/(9-1) 



(37) 



(38) 



As in the argument leading to Eq. (g) the fluctuations about the saddle point are canceled by the pre-exponential 
factor in Eq. ( pT| ) to leave power- law ('In-free') behavior 

The thermodynamic limit of the near-critical free energy, defined by Eq. (O) , now follows as 



where we have identified 



1 - a 



and (given Eq. 38) 



0+ 



(39) 



(40) 



(41) 



To establish the role of the remaining constants {r±) in Eq. ( |36| ) we consider the anomalous contribution to the 
free energy defined by 



F,(/3) ^ hm 1 F{f3, L) - lim 



Appealing to Eq. (^9|) , and recalling our sign convention [gO[ , we identify 



(42) 



(43) 



On the basis of rather general arguments po[ | we expect that away from a critical point the free energy anomaly is 
just minus the logarithm of the number of coexisting phases, so that 
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r+ = -Fa{p < /3c) = (44a) 
r_ = -Fa(/3>/?c)=ln2 (44b) 



In the critical finite-size limit we find from Eq. (|14|) 

F(/3c,i) = i''/o(/3c)-lni(0) (45) 
The critical value of the free energy anomaly, defining the Privman -Fisher constant Uq JT^, follows as 

Uo = Fa{Pc)^-\^Z{0) (46) 

These predictions are testable to varying degrees through both the energy-dependence of the energy pdf and the 
temperature-dependence of the associated free energy. 

Figure 6 shows the results for the ratio of the specific heat amplitudes that follow from Eq. ( ^ ) when the measured 
decay of the critical energy pdf (Figs. 4, 5) is matched to the prediction (^6|), in conjunction with Eq. (^7|). We 
can expect the predictions and observations to match up only in a window of x values. Clearly, x must be large 
enough to lie within the thermodynamic critical region; but it must also not be so large that the associated energy 
lies outside the domain of validity of the basic scaling ansatz (Eq. ^) . The size of this window should increase with 
increasing system size. The location of this window on the x-axis may also be expected to be different for the positive 
and negative x-branches -if, as seems reasonable, one regards the correlation lerwth ^ (rather than x ^ e — ec or 
y ^ f3 — (ic) as a measure of criticality. This is, indeed, the view we have adopted Q. Thus Fig. 6 shows the results 
for the 'effective' amplitude ratio, obtained by fitting over ranges of x values, with each pair of (positive and negative 
ranges) being centered on the same value of z = L/^, used as the abscissa |3^. On the basis of this data we make 
the assignment = 0.575(10) which is to be compared with A+/A„ — 0.523(9) in reference and 0.567(16) 

in reference 125(1. 



In Fig. 7 we show the results for InZ(y) that follow (cf Eqs. (pT|), (34), (p35|)) from the measured energy pdfs, using 
Eq. (35). The latter determines Suffix) only to within an additive constant which must be fixed by appeal to the 



predicted value of either r+ or r_. We have chosen the latter so that Eq. (44b) is satisfied, by fiat. The motivation 
for this choice is that it provides us with an inherently more reliable estimate of the parameter (which, unlike rj-, 
is not known a priori). Since —Uq — InZ(O) is closer to r- than to r^ the function \nZ[y) converges more quickly 
to its y > asymptote than it does to its y < asymptote. Fixing r_ (the intercept of the y > Q asymptote) thus 
tethers the value assigned to [/q more effectively than fixing r^. As with the amplitude ratio considered above, the 
value assigned to Uq depends upon the range of y- values used in the fit to the anticipated asymptotic form (Eq. ^7|) . 
Again we have chosen to characterize the temperature range utilized through the value of the ratio z = L/£^\ again 
we can expect the analysis to be trustworthy only if it is based upon data lying within the thermodynamic-critical 
window. Our data (Fig. 8) do not allow a systematic analysis of the approach to the desired limit; but they provide 
the basis for the assignment Uo — —0.57(2) [Q. The assignment of the uncertainty limit is subjective but, we think, 
conservative. We note the close correspondence with the assignment {Uq = —0.57) emerging from an earlier study 
p6[ , similar in concept, but utilizing the distribution of the order parameter. However our assignment differs (in what 
would seem to be a statistically significant fashion) from the result Uq = —0.625(5) obtained by Mon on the basis 
of altogether different techniques. 

IV. MICROCANONICAL AND CANONICAL SPECIFIC HEATS 

A. Generalities 

Thus far we have focused on the implications of the microcanonical entropy for observations made in the canonical 
ensemble. We now turn to consider their implications for observations made within ensembles that are (or are 
approximations to) microcanonical. 

We will assume (in keeping with eg [^|J^]) that the temperature of a microcanonical system should be identified 
from the relation 

^ j91nr2(e,L) 
/3''^(e, L) = L-'' (47) 

This identification is certainly required in the thermodynamic limit; but in the context of systems of finite size it is, 
it seems, a matter of convention ||37|. 
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It is illuminating to link this temperature with canonical observables. Appealing to Eq. ( ]15| ) we may write 

__,ainP(e|/3,L) 



de 



(48) 



where (notwithstanding appearances to the contrary) the rhs depends on e but not (3. This result shows that the 
equation prescribing the microcanonical temperature for a given energy is just the inverse of the equation prescribing 
the most probable energy for a given temperature: 



f3 = l3^'{e,L) 



6-(/3,L) 



(49) 



By comparison, within the canonical ensemble itself, 'the' energy for a given temperature is customarily identified 
with the canonical mean: 



(50) 



Eqs. ([49| ) and (^0|) make it immediately plain that the energy-temperature relationships associated with the two 
ensembles will coincide to the extent that the canonical energy distribution is gaussian (and thus has coincident 
mean e'^'^ and mode e'^'^). This correspondence is guaranteed in the thermodynamic limit, but not (in general) when 
finite-size effects are significant. 

The energy-temperature relationships are most usually probed through their derivatives, the associated specific 
heats. In the microcanonical case 



de 



\nP{e\l3,L) 



where, again, the /^-dependence of the rhs is illusory. 
In the canonical case (appeahng to Eq. (Esh) 



c^^(/3,L) 



de''%f3,L) 

a/3 



9e2 



= L'^e'^^\f3, L) 



(51) 



(52) 



Like the two 'caloric equations of state' (Eqs. (^^ and (^0|)) these two specific heats are guaranteed to agree in the 
thermodynamic limit; but they differ (in general) in the finite-size critical regime to which we now turn. 



B. Scaling forms 



First, we examine the asymptotic scaling regime where the the 'background' contribution to the specific heat can 
be neglected. We will consider the consequences of the corrections associated with the latter in the following section. 

In the scaling regime where the canonical energy pdf may be represented by its scaling form (Eqs. (^6|) and (|lT 
Eq. (E^) can be rewritten in terms of the energy and temperature scaling variables (Eqs. (P) and (|l^)) as 



dliiPix\y) _ dSeffix) 



dx 



dx 



(53) 



where in the last step we have exercised the right to set y = (the result is independent of y) and have made use of 
Eqs. (|l7| ) and (34). The microcanonical specific heat (Eq. ^) follows in scaling form 



with 



\nP{x\y = 0) 


-1 


d'^Seffix) 


dx"^ 




dx'^ 



c^'^ix) 



The scaling form of the canonical specific heat follows in a similar fashion, using Eqs. (|19| ) and (|5 

C-(/3,L)~i"/X'2~ce(y) 



(54a) 



(54b) 



(55a) 
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with (Eq. H) 

c-(y)=i(2)(y) (55b) 

The forms of both the scahng functions c'"^(a;) and c'^'^iy) can be determined from the scahng form for the micro- 
canonical entropy (Fig. 5) or, equivalently, the critical canonical energy pdf (Fig. 4), established in the preceding 
section. They are compared in Fig. 9. In the microcanonical case we have used Eq. ( p3| ) to identify the microcanonical 
temperature y = y^'^{x) to be associated with a given value of the energy variable x. 

In the thermodynamic limit realised at large values of | y | the two functions are, necessarily, consistent with one 
another, and approach the asymptotic behavior implied by Eq. (^9|). In the finite-size-critical (small | y |) regime, 
however, clear differences between the two scaling forms are apparent. In particular, the microcanonical maximum 
exceeds the canonical maximum by some 10%. One can show (Appendix C) that this -the fact that the microcanonical 
maximum is the larger one- follows necessarily if the scaling function S{x) is concave. 

The two scaling functions cross very close to the point [y — 0) identifying the bulk critical temperature. One can 
see this already in published microcanonical data ||]; a similar 'ensemble-independence' has also been noted in studies 
of the 'gaussian ensemble' [^8|. We have been unable to see any deep reason for this correspondence; but we do not 
discount the possibility that there is one. 

Though clearly visible, the differences between the two scaling functions are smaller than suggested by existing MC 
data HJ. The next section explains why. 



C. Beyond scaling: the role of the 'background' 

To understand the behavior observed in MC studies of microcanonical behavior, we must allow for the 'corrections 
to scaling ' which, in the canonical ensemble, are reflected simply in the existence of the additive negative background 
contribution to the heat capacity; their signature in the microcanonical ensemble is more subtle. 

The differences between the canonical and microcanonical results is effectively a strongly anharmonic effect: in a 
system of finite size, critical fluctuations sample a region of the entropy surface sufhciently large that the variation of 
its curvature becomes significant. We can expose the consequences analytically within an anharmonic perturbation 
theory in the cumulants of the energy pdf. The calculation is straightforward and we describe it in outline only. 

We appeal to the cumulant representation of the energy pdf at some (general) temperature: 



P{e\p,L) = ^ r dre-^exp 



.n=l 



n 



(56) 



We expand perturbatively to first order in the fourth cumulant and to second order in the third. We evaluate the 
second derivative of the logarithm of this function, which determines (cf Eq. (pl|)) the microcanonical specific heat 
associated with a given energy density. We evaluate this function at the modal energy e — e(/3, L) associated with the 
chosen temperature, prescribed by the (perturbative) solution of the microcanonical caloric equation of state (Eq. ^8[ ). 
The result is 



= (/3,L) 1 



;W(/3,L)e(^)(/3,L)-[e(3)(/3,L)]- 
2[e(2)(/3,L)]' 



(57) 



Eq. ^ shows that the cumulant correction terms displayed in this equation arc 0{L^'^) in the thermodynamic 
limit, confirming the equality of microcanonical and canonical predictions in this limit. To see what happens in 
the finitc-sizc-critical region we focus (for simplicity) on the temperature Pm for which the canonical specific heat is 
maximal, identified by the solution of 



.(3) 



(/3,L)=0: 



dc'"'{f3,L) 
dp 



(58) 



At this temperature Eq. ( p7| ) simplifies to 



c'^^(e„,L) =c^^(/3„,L). 



2e(2)(/3„,L) 



(59) 
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where em = ^{Pm,L), and we have used Eq. (p^. Now we appeal to the scaHng forms for the cumulants (Eq. |9|) 
and fold in the effects of the additional non-scaling contribution to the second cumulant (Eq. ^8|) to conclude that 



2[x(^){ym)-9{L)] 



(60) 



This result makes clear (albeit perturbatively) that, in the finite-size-limited regime, the temperature- independent 
additive 'background constant' in the canonical specific heat (manifested in the parameter g{L)) does not simply 
translate into an additive energy-independent background in its microcanonical counterpart. 

To expose the implications for the difference between canonical and microcanonical specific heats we introduce the 
dimensionless parameter: 



c="(/3m,i) 



Then 



RjL) 
R{oo) 



iVrn) 



i(2)(y„)-5(L) 



2[x(^){ym)-9{L)y 



COc 



(61) 



(62) 



where R{oo) is the scaling limit of R(L). 

The significance of the background constant cqc -in particular, its sign- is now apparent. The negative value of 
this constant results in an amplification of the difference between the microcanonical and canonical results (at /3m )j 
to a degree that diminishes with increasing system size. This is not simply the trivial effect that would arise from a 
uniform (downward) shift of both functions: Eq. ( |60|) shows that this is not what happens, as does the power of two 
on the rhs of Eq. (p3). 

It is not hard to track down the origins of this effect. The difference between the canonical and microcanonical 
specific heats is, we have noted, an anharmonic effect; in the present context the 'corrections to scaling' reduce (only) 
the second cumulant of the energy pdf and thus, in a relative sense, enhance the anharmonic (non-gaussian) character 
of the energy pdf, as one can see immediately from a comparison of Figures 3 and 4. 

The effect is significant. For L — 10 (as used in the simulations reported in ||]), estimating c'^'^{f3m, L) by c'^^{Pc, L) 
one can read off from Fig. 2 that R{L)/R{oo) ~ 4. The somewhat unexpected conclusion that the fractional difference 
between c'^^ and c'^^ at bulk criticality actually decreases for increasing L is consistent with some MC studies [B9[. 



V. CONCLUSIONS 



We review, briefly, the three principal strands of this work. 

First, we have broached the general question of the finite-size corrections to the density of states of a many-body 
system. The explicit proposal for the pre-exponential structure advanced in Eq. (^) is consistent with the pre- 
factor-free structure of the canonical partition function |10| and with the behavior of the simple models discussed in 
Appendix B. Given the growing interest in the behavior of mcsoscopically-sized systems, this proposal seems to merit 
some further study, with more rigor than we have attempted to offer here. 

Second, we have shown how one can fold out from the canonical energy fluctuation spectrum the principal corrections 
to scaling. The underlying behavior exhibits scale-invariance to a degree that seems remarkable, given the relative 
weakness of energy fluctuations. It is, we have seen, also largely consistent with established 3d-Ising critical properties. 

Third, we have provided a finite-size scaling theory of the microcanonical ensemble. This was the original motivation 
for this work- specifically, the suggestion |^ that the finite-size-smearing of critical behavior characteristic of the 
canonical ensemble is 'greatly reduced' within the microcanonical framework. Reference offers two pieces of 
supporting evidence for this contention, which merit final comment. 

Reference [ pf s uggests, firstly, that, in the vicinity of Cc, the microcanonical entropy (measured with the techniques 
described in ||40|) can be adequately represented by a form (Eq. 6 of reference |^) which allows for no finite-size 
corrections at all, and which corresponds essentially to the large x limit (Eq. |3^) of our scaling function. In fact the 
quality of the fit provided by this representation is rather poor. And we would expect it to be so. The measured 
microcanonical entropy evolves in a manifestly smooth way |4l| between the limiting thermodynamic forms appropriate 
above and below e^, Eq. 6 of reference H] is non-analytic at ec- Moreover, in analyzing data for the entropy and 
its derivatives, it is - we have seen - essential (on all systems practically accessible) to do justice to the corrections 
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associated with the background constant cqc- Even in the thermodynamic hmit the corrections allowed for in Eq. (8) 
of reference |^ do not do this. 

The second piece of supporting evidence offered in reference Q is a striking enhancement of the critical peak in 
the microcanonical specific heat, with respect to its canonical counterpart. As we have seen, this behavior is at least 
partly due to the effects associated with cqc; Figure 9 indicates that the underlying differences are rather less 
dramatic. 
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FIG. 1. The canonical mean of the energy density for the critical d=3 Ising model as a function of system size [Q. The 
statistical errors are an order of magnitude smaller than the system size. The points marked o are taken from reference [ p5[ . 
The relevant parameters have been assigned the values |2^]: a — 0.108, v = 0.63067 and /3c = 0.2216544. The arrow identifies 
the best-fit value for the intercept, prescribing the constant Ec (Eq. pl[ ). 
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FIG. 2. The canonical variance of the energy density for the critical d=3 Ising model as a function of system size [Q. The 
statistical errors are an order of magnitude smaller than the system size. The parameters are as specified in Fig. 1. The arrow 
identifies the best fit value for the intercept, prescribing the constant coc (Eq. pd). 
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FIG. 3. The canonical pdf (Eqs. |l| - 0) of the scaled dimensionless energy density x of the critical d=3 Ising model, for 
a range of system sizes. The scaled variable is defined in Eq. ^ with the choice (cf Fig. 1) tc = —0.9909. The scale factor Oe 
implicit in the scale of the x- variable is chosen such that (cf Eq. (0)) x = e — for L = 10. 
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FIG. 6. Estimates of the specific heat amplitude ratio, deduced from the decay of the energy pdf (for different system sizes) 
at 'large' (positive and negative) a-values. The estimates were determined by fitting to pairs of ranges of x values, with the 
ranges forming each pair being centered on a common value of L/^, which forms the abscissa. 
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FIG. 7. The function In Z{y) defined in Eq. (0), deduced from the microcanonical entropy (Fig. 5). The dimensionless 
variable y provides a scaled representation of the reduced (inverse) temperature (Eq. ^) . The str aigh t lines represent fits to 
the predicted asymptotic forms (Eq. The arrows identify the roles of the parameters r± (Eqs. 44a ,b) and Uo (Eq. 
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FIG. 8. Estimates of (Eq. 
the range defines the abscissa. 
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FI G. 9 . Comparison of the dimensionless microcanonical and canonical specific heat scaling functions, c*"^ (Eq. 54a) and 
(Eq. 55a), plotted as a function of the scaled dimensionless temperature (Eq ^). The light dashed line shows the power- 
behavior characterizing the thermodynamic (large \ y \) limit, extended back into the finite-size-limited region. 
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APPENDIX A: DEFINING A DENSITY OF STATES 



We discuss here, in general terms, the issues arising in defining a density of states function for a system in which 
the energy spectrum is discrete. The conventional argument ]43| makes the identification 

n{e,L)de^ (Al) 

E<E^<E+SE 

with the implicit assumption that the right hand side is proportional to SE (= L'^Se). This requires that: 

[CI] There exist many distinct levels s within the interval SE. 

[C2] The level degeneracy ilg is slowly varying over the interval E ^ E + 5E. 

The fractional variation of fig over the interval 5E can be estimated using Eq. (^7|); condition C2 can then be expressed 
in the form 

'^-}^^^5e = P^^{e,L)6E«l (A2) 
ae 

Taken together, conditions CI and C2 thus amount to the requirement that 

« -P^) ^^'^ 

where SEj characterizes the intrinsic discreteness of the energy spectrum. This condition is trivially satisfied in the 
classical limit (Appendix Bl considers one case explicitly). But there are obvious exceptions: in the Ising model 
(Appendix B2) Eq. ( |A3| ) is satisfied only at energies corresponding to microcanonical temperatures that are 'high' 
on the scale of the critical temperature. Or, to put it another way, ^Is is certainly not slowly vary ing over a range 
wide enough to embrace many system energy levels. We must now recognize, however, that Eq. (Al) (along with 
its implicit assumptions) does not faithfully reflect the conditi ons needed to legitimize the transition from discrete 
(Eq. la) to continuum (Eq. lb) representations. Instead of Eq. ( |Al| ) we require, rather, that we can consistently write 

n{e,L)e-'^^de^ ^ n^e'^^' (A4) 

E<E,<E+SE 

where (while retaining condition CI) we must replace condition C2 by 
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[C2A] rise is slowly varying over the interval E E + 5E, if thai interval lies in a range contributing significantly 
to the thermal properties at temperature 

The range contri buti ng significantly ... is centered on the saddle point E = L'^e (Eq. As a result, while condition 
C2 requires Eq. (|A2| ), condition C2A requires only that 

<P ln »(e, L) ^ ^1 ^ ^ 

where we have used Eq. (pi]). Thus, in place of Eq. ( |A3| ), wc need simply 

5Ei < [L'^c''^(e, L)] ^'^ ~ [e^^) {|3^'^ ^ L)j ^'^ (A6) 

where the last step uses Eq. (^), and /J^'^ = (3'^'^{e,L) (Eq. ^7|). This equation expresses more explicitly the implica- 
tions of condition C2A. A density of states function will exist in the operational sense (Eq. [lb| ) that it may be used to 
compute thermal properties at a given temperature as long as the canonical energy distribution (for that temperature) 
is broad on the scale of the intrinsic discreteness of the energy spectrum. 

APPENDIX B: DENSITY OF STATES OF SIMPLE MODELS 

Here we show that the general form for the density of states function proposed in Eq. (||) is consistent with exact 
results for two simple models. 

1. Quasi-continuous energy spectrum: harmonic lattice model 

Consider a system (a harmonic model of the vibrations of a crystal structure, for example) whose energy spectrum 
is that of N weakly-interacting harmonic oscillators, with associated frequencies Vj, j — I . . . N . Then 

N N 

E{{n}) = h'^UjVj = 

gives the energy of a microstate in which (for each j) mode j has quantum number Uj. We co nsider the classical 



(h 0) limit, in which the energy levels are quasicontinuous. In this case 5Ei ^ hvmini Eq- (A3) is satisfied, and we 



may proceed as in Eq. (Al) to write 



where 

DiX) 



1 if < A < NSe 
otherwise 



while e = E/N is the energy per oscillator. In the ft, — > limit the sums on Uj can be replaced by integrals on ej to 
give 

n{e,N) = e-^E^Q{N)I{N) 

where 

N ^1 



/(A) = f dx,<5(l-^a;,) 



and 
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Writing an integral representation of the (5-function wc find 



I{N) 



271 



dhe 



-ih 



Ah 



dhcos [{N - 2)h] 



sin h 



N 



which may be approximated using the asymptotic expansion for the F function 



(Bl) 



In analyzing the remaining (energy-independent) contribution we suppose that the frequency spectrum is that of a 
d = 1 system of particles, with a gap. Then 

InQ(iV) = -Nhih- Nqn 
where the sum may be written in the form 

N N-1 



QN 



where H{6) is periodic, and (invoking the assumed gap) H{0) is non-zero. It can be shown [ p5[ that the — > oo limit 
of this sum 



9o 



1 

2^ 



H{e)de 



has finite-size corrections that are exponentially small in N. 

Gathering these results together we conclude that the density of states has the form of Eq. (||) with the identifications 
N ^L'^ and 



1 



s(e) = In 

Appealing to Eq. (H) one can readily recover, as a check, the canonical partition function 



(B2) 



2. Discrete energy spectrum: Id Ising model 

Consider a d = 1 Ising model of N sites, with periodic boundary conditions. Choosing the ground state as the 
energy-zero, the energy density for a macrostate of M domain walls is e = Mei/N, where e/ is the domain wall energy. 
The number of microstates corresponding to macrostate M, N is 



2x m 



M,N 



(iV- M)!M! 



Appealing to the asymptotic form (Bl) once more we find that 



where x = M/N = e/ej. In this case SEj — e/ and Eq. (A3) is not in general satisfied. But, since Eq. (A6) is, we 
may still identify a density of states by 



n{e,N) = ^nM,N 
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which one may then readily recast in the form of Eq. (||) with the identifications N = L'^ and 

s{e) = -{I- x)\n{l- x) - x\nx (B3) 
Again, as a check, one can readily use this result to recover the canonical free energy density in the form 

/(/3,7V) = -^lnZ(/3,7V) = -\n{2coshK)+ K 

where 2K = Pej, and the last term reflects our choice of ground state energy. 

APPENDIX C: A BOUND ON THE CANONICAL SPECIFIC HEAT 

We outline here an argument establishing that the maximum of the value of the microcanonical specific heat 
provides an upper bound for the canonical specific heat, within the asymptotic scaling region. The argument assumes 
the concavity of the function Seffix); the concavity of SJx) is already presupposed in the formulation of Eq. (|^). 

We write the scaling function for the energy pdf (Eq. ^_7|) in the form: 

Pix\y) = Qix, y)G{x ~ x{y),~l/S':fj{x*)) (CI) 
where G{z, b) is a gaussian of zero mean, and variance 6; x{y) is the modal scaled energy, for a given y, the solution of 



dx 



and x* locates the maximum of the microcanonical specific heat, identified by the condition (Eq. 54b) 

~s':ff{x)<~s':fj{x*) 

The function Q{x,y) introduced in Eq. ( |Cl| ) is defined by: 

g(x,y) = Qoe^(^'^) 

where 

T{x, y) = -xy + Seff{x) - \ x - x{y)f 



(C2) 



(C3) 



(C4) 



(C4) 



while Qo is an x- independent constant, defined by normalization. From the assumed concavity of Seff{x) it is 
straightforward to show that, for an y giv en y, T{x, y) is concave in x, with a single maximum at a; = x{y). 
Now appealing to Eqs. (|Cl| ) and ( ^5b| ) we can write 



(C5) 



dxP{x\y) X — x'^^\y) 
< I dxP{x\y) [x - x{y)f 

dzQ{x{y) + z, y)G{z, -l/5^'^.^(x'^))z2 
dzQ{z,y)G(z,-llS'lff{x*))z^ 



(C6) 



where 



Qiz, y) = 77 [Q{x{y) + z,y) + Q{x{y) - z, y)] 



^T{x{v)+z,y) ^ ^T(x{y)-z.,v) 



(C7) 



26 



From the properties of the function T{x, y) it is straightforward to show that Q{z, y) has a single turning point (at 
z = 0), and that there exists some 2:0(2/) such that 



Q{z,y) 



> 1 if I z |< zo{y) 
< 1 if I z |> zo{y) 



(C8) 



Then, finally, appealing to Eqs. (C6) and (541:) 



2-(2/) - 5'-(x*) < I dz g(z,y)-l G(z, (^*))-^' 



< zliy) J dz [Q{z,y) - ij G(z, -l/S^'^/x*)) 

= 



(C9) 



where the last step exploits normalization conditions. 

It follows that the microcanonical specific heat maximum c^"^(a;*) provides an upper bound for the canonical specific 
heat. 
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